Critical coupling for chiral symmetry breaking in QCD motivated models 
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We determine the critical coupling constant above which dynamical chiral symmetry breaking occurs 
in a class of QCD motivated models where the gluon propagator has an enhanced infrared behavior. 
Using methods of bifurcation theory we find that the critical value of the coupling constant is always 
smaller than the one obtained for QCD. 



I. INTRODUCTION 

The hadronic physics at low energy is expected to be 
described by the infrared properties of Quantum Chro- 
modynamics (QCD), whose non-perturbative character 
in general forces us to make use of approximate models 
in order to understand the strong interaction effects at 
a small momentum scale. One of these models became 
known as Global Color Model (CGM) § - §], which is 
a quark-gluon quantum field theory that describes QCD 
for low energy processes. There are many recent calcu- 
lations exemplifying the remarkable success of this pro- 
cedure l^,!). It relates the hadronic properties to the 
Schwinger functions of quarks and gluons. Therefore, 
when comparing the theoretical calculations to some low 
energy data, e.g. pseudoscalars masses and decay con- 
stants or other chiral parameters, we are learning about 
the infrared behavior of the quark and gluon propaga- 
tors. In the near future, this semi-phenomenological tool 
may reveal to be even more successful than the relativis- 
tic quark model or the bag model. 

In order to have an idea of what is behind the GCM, we 
can recall that its action is obtained from the QCD gen- 
erating functional in the standard way with the mam 
difference being that in the functional generator of the 
connected gluon n-point functions we neglect the higher 
than two n-point functions, expecting that a phenomeno- 
logical description of the gluon propagator g^D^^{x — y) 
retains most of the information about the non-Abelian 
character of QCD. The effect of this approximation can 
only be measured in comparison with experiments, and in 
fact it does work quite well once we model appropriately 
the infrared behavior of the gluon propagator. 

Some of the gluon propagators used in the GCM have 
a quite enhanced infrared behavior. One example is the 
introduction of a delta function 5{k) as prescribed in 
Ref. 1^ , which is a confining propagator according to the 
criterion of absence of real fc^ poles for the quark propa- 
gators . Another ansatz for the gluon two-point func- 
tion in the infrared is g'^D{k^) — 37r^(x/A)^ exp (— fc^/A) 

, which was inspired by and approaches the 6 function 
ansatz of Ref. [|| for A ^ 0, where x ^^'^ A are ad- 
justable parameters. In detailed calculations of chiral 



parameters more complete ansatze for the gluon propa- 
gators than the above ones have been used, in general 
including the asymptotic behavior of the gluon propaga- 
tor as predicted by QCD. In Ref. ||] the following form 
for the gluon propagator in the Landau gauge was intro- 
duced 



D{k^) , 



(1) 



where, 

D{e) 



/c2[i-i-n(fc2)] 



4.71"^ rr?^ 5^ {k) 



1 _ p(-fcV[4m?I) 
P 



(2) 



with d = 12/(33 — 2n/), and Uf — i (considering only 
three quark flavors). The mass scale rrit determined in 
Ref. 1^ was interpreted as marking the transition be- 
tween the perturbative and nonperturbative domains. 
Another ansatz is 



g'Dik^ 



2 X 



exp (-fc^) + 



fc2 



where 



a,(fc2) 



47r2rfg2 



(3) 



(4) 



J-{k'^)is a function chosen differently in the papers of 
Ref. 1^, as(fc2) describes the QCD running coupling con- 
stant where r is a parameter adjusted phenomenologi- 
cally. Notice that in the above expressions the momenta 
are in Euclidean space. Of course, there are still other 
variations of these ansatze |10| and attempts to explain 
the enhanced behavior in the infrared [pT| . 

In this work we study the bifurcation of the quark self- 
energy within the context of the GCM, i.e. we determine 
the critical coupling constant of the truncated Schwinger- 
Dyson equation for the quark propagator above which 
dynamical chiral symmetry breaking occurs using the 
gluon propagators discussed in the previous paragraph. 
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It is known from analytical and numerical studies of the 
Schwinger-Dyson equations that dynamical chiral sym- 
metry breaking takes place in QCD when the coupling 
constant (as) is of 0(1) fl^-IT^. These studies have also 
been performed with the use of effective potentials jl^ ] 
and corroborated by lattice simulations [Q. Therefore, 
it is natural to ask what is the critical coupling in the 
GCM with the phenomenological gluon propagators pro- 
posed in the literature. 

We will apply the standard techniques of bifurcation 
theory as used by Atkinson and collaborators [ p^ , p7| pj|] 
and verify that for the gluon propagator given by a delta 
function the chiral symmetry is always broken no mat- 
ter what "perturbative" propagator is added to the 6. 
For a propagator of the form given by Eq. we obtain a 
lower bound on the critical coupling smaller than the one 
obtained for QCD and evaluate the smallest char- 
acteristic number confirming the result obtained with 
the S function propagator in the limit that the gaussian 
width goes to zero. In Section II we discuss chiral sym- 
metry breaking for an infrared gluon propagator similar 
to Eq.(|^) using a bifurcation analysis of the Schwinger- 
Dyson equations. In Section III we analyse the infrared 
gluon propagator given by a gaussian form using a differ- 
ent technique. The last section contains our conclusions. 



II. THE CRITICAL COUPLING FOR A 
CONFINING PROPAGATOR 

The Schwinger-Dyson equation for a massless quark 
propagator in Minkowski space can be written in the form 



S-\p) ^^-ig' 



2 



2 ' 



(5) 



where the gluon propagator in the Landau gauge, which 
will be used throughout the paper, is given by 



cab 



fc2 



(6) 



In the case that the gluon propagator is given by the 
confining delta function |q| 



9)^1' 



(36^ (k), 



(7) 



where /3 is an adjustable dimensional parameter. Note 
that the term k^ki,5{k)/k'^ is undefined from the point of 
view of generalized functions. Of course, we could work 
with the Feynman gauge where this unwanted behavior 
is softned. Fortunately the 5 function gives an integrable 
singularity and only because of this peculiarity we have 
not a strong pathological behavior. With Eq.(|^) it is 
quite easy to verify that the quark self-energy has a non- 
trivial solution for any positive value of g^, in the rainbow 



approximation (r'^(g,p) =7''). With the quark propa- 
gator given by 



S-\p)^A{p^)^~B{p% 



(8) 



and the gluon propagator of Eq.(Q) it follows from Eq.(||), 
in Euclidean space, the set of coupled integral equations: 



W) - l]p^ 



(2^ 



S\p-q) 



Aiq^ 



p-q + 2 



q^A^{q^)+B'^{q^) 

1- {p-q){p-q) -p 



d^q 
{2ny 



■5\p-q) 



[p - 

B(q^) 



q^A'^{q^) + B^{q^ 



The solutions of the nonlinear coupled integral equa- 
tions were determined in Ref. 0J^, showing that v4(p^) 
is a constant. Without loss of generality we can set 
A{p'^) = 1 and verify the behavior of the self-energy apart 
from a normalization constant. In this case we obtain the 
following equation for BijP') 



Bip') 



f^, B{q^)5\q-p) 



{2nY 

whose integration leads to 



52 (q2) 



B{p') 



B{p^) 



(27r)4p2+52(p2) 

Equation (O) has the solution 



(10) 



(11) 



Bip") = 



p2 < 







ri2 ^ 4£f£ 



for any non-zero and positive value of the coupling con- 
stant. If we consider only this infrared behavior of the 
gluon propagator we would conclude that this model is 
inconsistent with previous work on QCD, which predicts 
a critical value of the coupling constant for the onset of 
chiral symmetry breaking. However, the ultraviolet part 
of the gluon propagator cannot be neglected and the use 
of the following propagator is frequent |^J^ 



f3S\k) 



k^ky 

fc2 



1 

A?2' 



(12) 



The following comments are in order: a) In some papers 
the coupling constant is absorbed in the parameter 
(} ||]. Here we assume that g"^ can always be factorized, 
and the critical behavior is going to be related to it. b) 
We assume a constant g'^ . In the next section it will be 
shown that our result does not change with the inclu- 
sion of the running coupling in the ultraviolet part of the 
gluon propagator, c) This propagator has been used in 



(9) 
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the literature with some minor differences; these varia- 
tions are not important for our result as will become clear 
in the next section. 

When the gluon propagator is given by Eq.(^) we 
need a thorough analysis of the quark self-energy to find 
for which value of the coupling constant the quark self- 
energy admits a nontrivial solution. The substitution of 
Eq.Q into Eq.(||) yields 



B{x) 
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47r2 



kB{x) 



^ dy yB{y) 

+ BHy)J' 



(13) 



where x =p^ ,y = q^, x^ax = max(a; y) and k= -^.^e 
will be looking for solutions of Eq. ( [l3| ) only in the interval 
[/i^,A^], although the 5 function has been integrated in 
the full range of momenta. A rigorous handling of this 
integration does not change the conclusions about the 
bifurcation point in the limit of very small (large) infrared 
(ultraviolet) cutoff. 

It is important to verify that Eq. ( p^ ) is a very peculiar 
one, in the sense that for very small values of 5^, keeping 
the product g^K fixed, we do have the solution of the 
nonlinear equation, which, at low momenta, approaches 




(14) 



^0 . K f ixed 



To find nontrivial small solutions of the nonlinear equa- 
tion (^3|) we study the linearized equation, i.e. the 
functional derivative of Eq.(|l3|) evaluated at B{x) = 
P,|lM9|. Writing 5B = /, we obtain 



dy 



(15) 



Note that we are not considering the equation for v4(p^) 
because this equation is of second order in 5B ||l^ . The 
existence of a solution for Eq. (|l^) is a sufficient condition 
for the onset of chiral symmetry breaking |l^ . Defining 
a = 5^/47r^ = as/vr, Eq.(jl^) is equivalent to the differ- 
ential equation 



x{x — Q!k)/ {x) + 2xf [x] + af{x) = 0, 



(16) 



together with the ultraviolet boundary condition [{x — 
Q!K)/(a;)]'|^_,A2 = and another condition valid in the 
infrared region. 

The choice of the infrared boundary condition is a cru- 
cial one. Note that if we neglect the second term of the 
right-hand side of Eq.(13) we only obtain a trivial result, 
no matter what is the boundary condition coming from 
Eq.(13). However, we know that the nonlinear equation 
always has a nontrivial solution for the 5 function prop- 
agator. The most suitable infrared boundary condition 
is the one that comes from the nonlinear integral equa- 
tion (12), since when applied to the linear equation the 



result must be consistent with the known solution of the 
nonlinear equation for small (see Eq.(|l4[)). Therefore, 
our infrared boundary condition is given by 



(17) 



Equation ( [iq ) can be put in the form of a hypergeomet- 
ric equation performing a shift x — > z -I- an and defining 
y = —z/an leading to 

2/(1 - y)f {-any + an) + 2(1 - y)f {-any + an) 

+Q!/(— Q!Ky + q;k) = 0, (18) 



The solution of Eq. ( |lq ) that obeys the infrared bound- 
ary condition is 



f{,x)=AF (a,&;2;l 



X 

an 



where 



A = Va^r(3/2 - CT)r(3/2 + a), 



(19) 



(20) 



is a normalization constant, a is defined as cr = (j—a)^^'^, 
and 



1 

2 ' 

1 

2 



1 



a. 



(21) 



To apply the ultraviolet boundary condition we recall the 
following relation involving hypergeometric functions 

-^[z^-^Fia, b; c; z)] = (c - l)z''-^F{a, 5; c - 1; z), (22) 
dz 

which lead us to the analysis of the zeros of the following 
equation 



[{x - aK)/(a;)]'U^A2 



- anAF 1/2 + CT, 1/2 - cr; 1, 1 



(23) 



+A2 



We need to consider the solution of Eq.([2^) in three 
different regions of the parameter a, namely, < a < ;|, 
a = J and a > j, and study the asymptotic behavior of 
f{x) in each one of these cases. When a < j the relation 



F(a, b; c; z) - 
Tic)Tib~a) 

r(6)r(c-a) 
, r(c)r(a-6) 
r(a)r(c-6) 

can be used together with 

F(a,6;c;0) = 1 

to give 



(-z)-"F(a, l + a~c;l + a-b; z'^) 
{-z)-''F{l + b-c,b;l + b-a;z-'^), (24) 



(25) 
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[{x - aK)/(x)]'|^^A2 = 

^ X-l/2+<T 

A' ( — ~ 1 I +B' 

an 



X 

an 



-1/2-0- 



(26) 



where A' and B' depend on a and an. 

In this case cr is a real and positive number smaller than 



For large values of x we see that Eq.(26) condition is 



satisfied. Therefore, ( pj| ) is a solution when < a < j. 

Note that when a = ^ we have (7 = 0, the constants 
a and b are identical, implying that ( p^ ) cannot be used 
in this case. However, we can perform the limit a; — > 
already in the differential equation ( p^ ) and study its 
behavior for this particular value of a. In the asymptotic 
region (UV) we obtain 



J{x)uv =x~^''^{C + D\nx). 



(27) 



We easily see that the ultraviolet boundary condition is 
satisfied for large x and (HOh is also a solution when a = 



When a > T in the parameters a and h of Eq.(19) we 



make the substitution a ip with p defined as (a — -j) 2 . 
The asymptotic behavior in this case is obtained with the 
help of (ll) and ^ leading to 

[{X - aK)/(.T)]'|j.^A2 = 

"r(2ip)r(3/2 - cr)r(3/2 + a) 



X 

XI 1 

an 



(28) 



which is valid for x ^ an. Eq.(|2^) has an infinite set of 
zeros located at a; = an^Xn + 1) for integer n, which can 
be determined with the procedure of the second paper of 
Ref. 0. 

The real part of Eq. (ESh can be written as 



exp < In 



^(2^p)^(3/2-g)^(3/2 + g) -1- 



(29) 



as Inz — ln|z| + i(argz ± 2m:), where z is a complex 
number and n = 0, 1, 2, . . ., the zeros of this function will 
occur for 



ln(a::„ 



2n7r — arg 



/r(2ip)r(3/2 - CT)r(3/2 + a) 



r2(i + ip) 



(30) 



with n = 1, 2, . . .. This relation was obtained for x S> an, 
and n — was excluded in order to obtain only positive 



values in Eq.(30). Therefore, also for a > i we do have 



a nontrivial solution. 

In summary, when the gluon propagator is modeled by 
a delta function plus a "perturbative" propagator 
we verified that the chiral symmetry is broken for any 
positive value of the coupling constant. One may ask 



if this result changes if we modify the non-perturbative 
as well as the perturbative propagator, as in the case of 
Eq.(^, or with the inclusion of the running coupling or 
the effect of a dynamical gluon mass [ [l3|j2(]f| . We will 
show in the next section that this is not the case. 



III. AN INFRARED MODEL IN THE GAUSSIAN 
FORM 

As discussed in the introduction a different ansatz for 
the gluon propagator, exemplified by Eq.(|^), has fre- 
quently been used. Its non-perturbative infrared behav- 
ior is given by the first term of Eq. (||) 



ra6 



(31) 



where the parameter x controls the intensity of the inter- 
action, and A gives the gaussian width. It is important 
to note that in the limit A — + we recover the infrared 
behavior of the propagator discussed in the previous sec- 
tion. 

In this section we will show that the quark self-energy 
calculated with the propagator of Eq.(^l]) has a well 
defined bifurcation point. Afterwards, we verify that 
adding a perturbative tail to this propagator results in a 
smaller critical coupling constant than the one obtained 
only with the perturbative part. Finally, in the limit 
A 0, we recover the result of the previous section 
showing that it is independent of the perturbative part 
that is added to the delta function. 

To find nontrivial small solutions of the nonlinear 
Schwinger-Dyson equation for the quark self-energy we 
can study the linearized equation. The substitution of 
Eq.(^ll) into Eq.(||) gives the following linear integral 
equation for f{x) (recalhng that f{x) — 5B{x)) 



fix) 



3 2X' 
4^ A^ 



dy/(y)exp — 
^ 



dyfiy) exp 



A _ 
(32) 



where we considered A{p'^) = 1, since, according to bi- 
furcation theory, the equation for (A(p^) — 1) is of higher 
order in the functional derivative of B{p^), and the criti- 
cal coupling constant is determined only through Eq. (|2|) . 
To obtain Eq. ( ^ ) we performed in the gluon propagator 
the so called angle approximation, which is given by 

D{{p - qf) « 0(p2 _ q^)D{p^) + e{q^ - p')D{q'). (33) 



The Eq. (^2|) is a homogeneous Fredholm equation with 
the kernel 



{x, y) = exp - y) + exp (-■^) ^{y - x). 

(34) 



The norm of (p3) is easily calculated 



4 



Kj 



dx 



dyK]{x,y) 



A2 
2 ' 



(35) 



and we find nontrivial solutions of B{x) for on a 
point set whose smallest positive point satisfies [ [17| 



9'> 



4A2 
3^ 



iff 



A 



(36) 



Eq.(p6[) gives a lower bound for the critical point g^. 
However, we can do better than this. Using the method 
of traces we can show that the approximate critical value 
is indeed of the order of the smaller value of the bound 
given by Eg. (^ . 

The following approximate formula holds true for the 
smallest characteristic number [pi| 



(37) 




where for a symmetric kernel 



A 



2m 



Kl,{x,y)dydx, 



with m running over 1, 2. K2{x, y) is given by 



K2{x,y)^ / Ki{x,z)Ki{z,y)dz, 

JQ 



(38) 



(39) 



and, in the case of the infrared part of the propaga- 
tor, Kiix,y) = {ix^/4,/S?)g'^Ki{x,y). The calculation 
of Eq. (|37|) entails 



' 9xV32A2 
297xV4096A4 



A 

1.97—. 



(40) 



It is known that this method overestimates g^. In this 
case we verify that the smallest value of our lower bound 
Eq.(^6|) is a good approximation for the critical point. 
The value of the critical coupling can be obtained sub- 
stituting the phenomenological values of x ^'^d. A into 
Eq. ( ^6| ) . It is obvious that in the limit A ^ we recover 
the result of the previous section, i.e. with the gluon 
propagator given hy & 5 function the chiral symmetry is 
broken for any coupling constant g^ > Q. 

Let us now consider the case where we add to the gaus- 
sian form of the infrared propagator a perturbative con- 
tribution. If this new contribution is of the form given 
by (H) or (H), it is not difficult to verify the existence of 
a critical coupling constant for the onset of chiral sym- 
metry breaking. The reason for this is that the ultravi- 
olet behavior of the propagator is softened by the fac- 
tors 1 — exp(— A:^/[4mt]) and T{k^), producing an effect 
equivalent to a dynamical gluon mass pO| for which it 
was shown the existence of a bifurcation point Jl3| , p0t . 
Therefore, we proceed to the most intriguing case where 
we add to the ansatz for the infrared gluon propagator 



the contribution of QCD with massive gluons and the ef- 
fect of the running coupling constant, i.e. a propagator 
proportional to as(A:^)/(fc^ + m^). 

If we denominate by Ku the kernel related to this ul- 
traviolet propagator, we assume that we do have a non- 
trivial solution for 



9'> 



1 



(41) 



A lower bound for the bifurcation point in the case that 
we consider the sum of propagators will be given by 



3x^ 
4A2 



Ki + Ku |1> 1, 



(42) 



where || Kj 
inequality 



is given by Eq. (|35|) . With the triangle 



II Ka + Kb \\<\\ Ka II + II Kb 
and making use of Eq. ( p5[ ) we obtain 

,2 



9 



2/3V2r 

I 8 A 



u 



> 1. 



(43) 



(44) 



The solutions that break the chiral symmetry appear 
only for values of the coupling constant obeying Eq. (|4|) . 
Note that as long as we can factor out g^ from both 
propagators, and as long as the kernel of the perturbative 
tail is bounded, we obtain a condition for values of g^ 
above which the self-energy bifurcates. The lower bound 
for the critical value 



9c 



> 



1 



II Ku 



(45) 



approaches zero as we take the limit A — > 0, and is com- 
patible with the result of the previous section. 

Considering the infrared and ultraviolet contributions 
to the propagator, we apply again the method of traces to 
determine the critical coupling constant through Eqs.(37 
and (p8|). In this case the kernel Ki{x,y) in Eq.(| 
given by 



Ki{x,y) ^ H{x)d{x ~ y) + H{y)9{y ^ x), 



where 



H{x) 



4A2 



exp 



A 



ml ln(^ 



(46) 



(47) 



Note that we have already factored out from the kernel 
the coupling g^. To determine Ai we need K2{x, y) which 
is equal to 

K2{x,y) ^ H{x)H{y) I dz + H{x) dzH{z) 

Jo Jy 

dzH^{z). (48) 



5 



The critical coupling constant was determined numer- 
ically for four quark flavors, = 2.0 GeV^, rrig = 600 
MeV, A = 300 MeV, and we assumed r = 16, where all 
these values are phenomenologically acceptable, and, in 
particular, the value of r is compatible with the infrared 
behavior of the running coupling constant in a theory 
admitting dynamically generated gluon masses [ p2[ (in 
these papers we have r « 4(mg/A^)). The critical cou- 
pling is shown in Fig.(l) as a function of the parameter 
A. Note that as we decrease A the value of the critical 
coupling goes to zero, confirming the result of the pre- 
vious section. In this region also the integrals appearing 
in Eq.(^) become more problematic because the gaus- 
sian becomes very peaked. Changes in the form of the 
ultraviolet propagator barely affect this result. Finally, 
we confirm that for this class of infrared propagators the 
chiral phase transition happens at a smaller value than 
the one obtained with perturbative QCD. 



behavior of the gluon propagator implies that the chiral 
symmetry is always broken for any value of the coupling 
constant. This is an interesting result if we remember 
that this propagator is considered to describe confined 
quarks Q. However, this is certainly in contrast with 
what is known to happen in QCD. 

With a model for the gluon propagator inspired by 
the delta function but softer at the origin (in a gaussian 
form), we verified, using the method of traces, that the 
critical coupling for the onset of chiral symmetry break- 
ing is lower than the one expected for perturbative QCD, 
and recovered in a particular limit the result obtained 
with the model involving a delta function for the infrared 
gluon propagator. Therefore, the critical point for chi- 
ral symmetry breaking may distinguish among different 
QCD motivated models, and in the cases we have stud- 
ied the critical coupling is always smaller than the one 
known for QCD MK^. 
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FIG. 1. Critical coupling constant, g^, as a function of the 
parameter A, considering the infrared and ultraviolet contri- 
butions to the gluon propagator. We see that decreasing the 
value of A the critical coupling goes to zero. The evaluation 
was perfomed for Uf = 4, = 2.0 GeV^, rUg = 600 MeV, 
A = 300 MeV and r = 16. 



IV. CONCLUSIONS 

An ansatz for the infrared gluon propagator in the form 
of a S*{q) function (or a variation thereof), frequently 
used in applications of the Global Color Model, gives 
an excelent description of hadronic properties connected 
with chiral symmetry breaking. In this work we studied 
the bifurcation of the quark self-energy within this model. 
The idea was to compare GCM to QCD with just the 
perturbative gluon propagator, where it is known that 
above a certain critical value of the coupling constant the 
chiral symmetry is broken JT3|-p^. We verified that the 
introduction of a delta function to describe the infrared 
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